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5-Meson Gateways to Missing Charmonium Levels 
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We outline a coherent strategy for exploring the four remaining narrow charmonium states 
[r;c(2^So), hc{\^Vi), ?7c2(l^D2), and V'2(1''D2)] expected to lie below charm threshold. Produced 
in B-meson decays, these levels should be identifiable now via striking radiative transitions among 
charmonium levels and in exclusive final states of kaons and pions. Their production and decay 
rates will provide much needed new tests for theoretical descriptions of heavy quarkoriia. 

PACS numbers: 14.40.Nd,14.40.Gx,13.25.Gv 



Every new spectroscopy carries the potential to il- 
luminate older spectroscopies in novel — and sometimes 
unexpected — ways. The decays of charmed mesons have 
offered new paths to the study of mesons — especially 
excited states — composed of light quarks, and decays 
of charmonium (cc) states have provided new access 
to glueballs — hadrons composed largely of gluons. The 
study of semileptonic r decays has considerably enriched 
our knowledge of ai properties. i3-meson decays into 
charmonium states are an indispensable tool for the ex- 
ploration of CP violation [1, 2]. They can also serve as 
gateways to the discovery of hitherto unobserved charmo- 
nium states. The properties of these missing states can 
illuminate the interquark force and reveal effects that lie 
outside the simple quarkonium potential framework, in- 
cluding the influence of virtual decay channels. 

Detailed knowledge of the cc spectrum first derived 
from the study of e"'"e~ annihilations, which explored the 
^Si and ^Di J^^ = 1 states and the ^Pj or ^Sq states 
connected to them by El or Ml radiative transitions. 
Though incisive, the e+e~ annihilation channel has lim- 
itations. Twenty years have passed without a confirma- 
tion of the Crystal Ball claim of the 2^80 77;.(3594±5) [3], 
and the complementary technique of charmonium forma- 
tion inpp annihilations does not support the 7y^(3594) ob- 
servation [4]. The sighting of the l^Pi level near the j 
center of gravity in pp — » /ic(3526) 7r° J/0 reported by 
Fermilab Experiment E760 [5] needs confirmation. The 
'7c2(l^D2) and V'2(1'^D2) have also proved elusive, with 
only an unsubstantiated claim in the literature for a 2 
state, V'(3836± 13), in vr^A^ J/^vt+tt- + anything [6]. 

Help is on the way. The CDF experiment's obser- 
vation of B decays to known charmonium levels shed 
light on prompt and secondary production of (cc) states, 
sharpening the puzzle of the production mechanism [7] . 
The e+e- ^ T(4S) ^ BB experiments CLEO [8], 
BaBar [9], and Belle [10], report increasingly detailed 
observations of the established quarkonium states in B 
decay. And now Belle has reported the observation of 



the 2^So (77^) level at a new mass (3654 MeV) in exclu- 
sive B K K sK~ -K^ decays [11]. 

In this paper, we give a template for the spectrum of 
charmonium based on a Coulomb + linear potential, and 
we present estimates of the principal decay rates for the 
unobserved states to make quantitative the expectation 
that four narrow states remain to be studied. We ar- 
gue that we may expect ample production of the missing 
charmonium states in i?-meson decays, and we suggest 
experimental strategies for detecting the missing levels. 
As with the discoveries of the P-states Xd and rjc of 
charmonium, radiative transitions will be of central im- 
portance. Finally, we comment on what we will learn by 
studying the masses and properties of the missing levels. 

The importance of radiative decays to the discovery of 
charmonium levels, including the D-wave states, has been 
appreciated since the earliest days of charmonium spec- 
troscopy [12]. Recently, Ko, Lee and Song [13] discussed 
the observation of the narrow D states by photonic and 
pionic transitions, and Suzuki [14] emphasized that the 
cascade decay B hcK^*^ — > ^rjcK'^*^ offers a promising 
technique to look for he. 

The spectrum of charmonium. To estimate the posi- 
tions of the missing charmonium levels, we have adjusted 
parameters of the classic Cornell potential [15, 16], 



(1) 



1.84 GeV, K = 0.61, 



a 



2.38 GeV" 



to reproduce the observed centers of gravity of the IS and 
IP states [38]. No one has produced a satisfactory an- 
alytic (or potential-based) explanation of the spin split- 
tings of the l^Pj levels. Moreover, the 2S, ID, and 2P 
levels are certain to be influenced appreciably by coupling 
to decay channels. Accordingly, we will not estimate the 
spin splittings of those levels beyond offering the expec- 
tation that they will be small. Our expectations for the 
charmonium spectrum are summarized in Table I. 

The fiPi he and the 2^80 v'c of course lie below DD 
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TABLE I: cc spectrum in the Coulomb + linear potential (1). 



State Mass (MeV) Remarks 



I'So 
l^Si 


Vc \ 


3067. < 


I 2979.8 ± 1.8 
1^ 3096.87 ± 0.04 


l^Pi 


Tic 


3526. 




l^Po 






C 3415.0 ± 0.8 


l^Pi 


Xcl > c.o.g. 


3526." < 


3510.51 ±0.12 


l^Pz 


Xc2 J 




[ 3556.18 ±0.13 




^ C.o.g. 


3678. < 


f 3654. ± 10. 
3685.96 ±0.09 


1^D2 




3815. 


-/^ DD (parity) 


l^Di 


V ] 




' 3769.9 ± 2.5 




1p2 > C.o.g. 


3815. < 


DD (parity) 




V-s J 




DD 


2P 




3968. 




IF 




4054. 




3S 




4118. 








3729.0 


t;ViTP<^linlrl 




D+D- 


3738.6 


threshold 


D°D*° or D*°D° 


3871.2 


threshold 






3879.3 


threshold 




D+D7 


3973.2 


threshold 






4013.4 


threshold 




D*+D*- 


4020.0 


threshold 


D+Dr or D*+D- 


4099.0 


threshold 




Dt+D*- 


4224.8 


threshold 



"Input values. 



threshold, and so wiU bo typically narrow charmonium 
states. In the absence of strong influence from the cou- 
pling to decay channels, the 2^Pj x'c ^'^d 2^Pi h'^ states 
should lie well above the DD and D* D thresholds, and so 
should have uninhibited strong decays. As has long been 
known, the J^'^ = 2"+ l^Dj 7?c2 and J^^ = 2" l^Dj 
'\\}2 states constitute an important special case: they lie 
between the DD and D*D thresholds, but are forbidden 
(because of their unnatural parity) to decay into DD. 
It is therefore plausible that they will appear as narrow 
levels, and we now quantify this suspicion. 

Properties of the missing levels. To estimate the de- 
cay rates, we shall use the established values for the 
r]c, J/tp,Xc,'>P' , and '0(3770) states, adopt the Belle value 
for M^/ , set Mh^ = 3526 MeV, and choose M^^^ = 
= 3815 MeV. We estimate the rates for hadronic 
and radiative decays in turn. 

Among hadronic decays, we consider transitions (tttt 
emission) and annihilations. To estimate the tt+tt" ± 
tt'^tt'^ transition rates, we use the standard multipole ex- 
pansion of the color gauge field [17, 18, 19] to express the 
El-El transition rates through the Wigner-Eckart theo- 
rem given in Eqn. (3.5) of Ref. [20], with experimental 
inputs given in Table X of that paper. The results are 
shown in Table II. For present purposes, the essential 
lesson is that we do not expect the tttt transition rates to 
be large for the missing levels of charmonium. 



TABLE II: Hadronic decay widths of charmonium states. 



cc state 


Decay 


Partial Width 


I'So 


Vc gg 


17.4 ± 2.8 MeV [21] 


l^Si 


ggg 


52.8 ± 5 keV [22] 


l^Pi 


he ggg 


720 ± 320 kcV 


l^Po 


Xco gg 


14.3 ±3.6 MeV*" 


l^Pi 


Xcl ggg 


0.64 ±0.10 MeV'' 


1='P2 


Xc2 gg 


1.71 ±0.21 MeV' 


2^80 


»?c gg 


8.3 ± 1.3 MeV= 




7]'c — > nnrjc 


160 keV* 


2^81 


i>' -» ggg 


23 ± 2.6 keV [22] 




tp' — > ■KlvJ/lp 


152 ± 17 keV [22] 




tp' rjj/ijj 


6.1 ± 1.1 keV [22] 


l^Da 


Vc2 gg 


110 keW^ 




??c2 mrrjc 


45 keV* 


l^Di 


V' ggg 


216 keV-'' 




ip — » ■wivJ/ip 


43 ± 15 keV 


l^Da 


i>2 ggg 


36 keV-'' 




tp2 — > mTJ/4> 


45 keV* 


l^Dg 


ip3 -> ggg 


102 keyf 




1p3 — > -WK-J/lp 


« 45 keV* 



"Computed from ^Pj rates using formalism of [23]; also see [24]. 

''Compilation of data analyzed by Maltoni, Ref. [23]. 

"^Scaled from T{ric — » gg). 

•^Computed using Eqn. (3.5) of Ref. [20]. 

"^Computed using Eqn. (3). 

■'^Computed using Eqn. (2). 

sProm rates compiled in Table X of Ref. [20]. 

For the annihilations into two or three gluons, we use 
the standard (lowest-order) perturbative QCD formu- 
las [25] to scale from available measurements for related 
states. This is a straightforward exercise for the S-wave 
levels. We use Maltoni's analysis [23] of the ^Pj annihi- 
lation rates to estimate the rate for he ggg. The rates 
for annihilations of the "^Dj states into three gluons (via 
color-singlet operators) are given by [26, 27] 

r(3D^ - ggg) = ^Cj ^-^^ lnAm,{r) , (2) 

where rI^} = d^R„i{r)/dr'^\^^^^, (r) = J^rrulf,{r), and 
Cj = 1, 4 for J = 3, 2, 1. A complete analysis (includ- 
ing color-octet operators as well) has too many unknowns 
to be of use [39]. The strengths of the J = 3,2, 1 anni- 
hilations are more generally proportional to Cj, even if 
color-octet operators dominate [28]. The two-gluon an- 
nihilation rate of the state is given by [29] 

3 m° 

Our estimates for the annihilation rates are collected in 
Table II. The expectation for T{ri[. gg) is to be com- 
pared with the Belle value of 15 ± 24 (stat) MeV [11]. 

The most prominent radiative decays of charmonium 
states are the El transitions, for which the rate [29, 30] 
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is given by 

r(n 



El 



-/ 2s+l«/ 



(2 J' + l)k^ 



(4) 



where Cc = ^ is the charm-quark charge, fc is the photon 
energy, the El transition matrix element is Sne-.n'e' = 
f f^runeir)un'f,ir) [f .7o(f ) - )] +0(fc/m,), n = 
n — £ is the radial quantum number, and {:::} is a 6-j 
symbol. For Ml transitions, the rate is given by 



/2s' + 1/, 



3m2 



(5) 

where Mnt-.n't = Unt{r)un' i{r)jo{^). 

The calculated rates for the prominent transitions 
among charmonium states are shown in Table III. Val- 
ues enclosed in parentheses have been corrected for the 
effects of coupling to decay channels, following the proce- 
dure developed in [16]. The calculated values reproduce 
the patterns exhibited by measurements, and are in good 
agreement with other calculations in the literature [31]. 
Wc expect them to provide reasonable guidance to the 
radiative decay rates of the missing charmonium levels. 

Integrating all the calculated rates, we note that the 
radiative decays should be prominent, with branching 
fractions B{hc ?7c7) « |, B{r]c2 hc'j) « |, and 
B{tp2 Xci,27) ~ |i of which B{tjj2 Xdl) ~ §• 

Charmonium production in B decays. Expectations 
for the fractions of S-meson decays leading to charmo- 
nium production are presented in Table IV. To estimate 
the B ^Sq production rates, wc appeal to the sug- 
gestion [35] that the ratio of spin-singlct to spin-triplet 
decay rates is relatively insensitive to poorly calculated 



-X) 



matrix elements, T{B n^Si+X)/T(B 
l + 8ml/ml 1.5. The inclusive production of IP states 
in B decays can be expressed [40] in terms of color-singlet 
and color-octet contributions as [36]: 



r(6 ^hc + x)/r{b + x) 

T{h^XcO+X)/T{h^£-Ui+X) 
T{h ^ xci + X)/V{b ^ rue + X) 
T{h^Xc2+X)/T{h^rue + X) 



U.lHs 
3.2H8 

UAHi + 9.3H8 
15.3#8 (6) 



Using the measured rates for inclusive Xci ajid Xc2 pro- 
duction summarized in Table IV we extract Hg, = (8.95 ± 
3.79) X 10-5 and Hi = (2.18 ± 0.31) x 10-^ which de- 
termine the inclusive branching fractions for XcO and he- 
No measurements exist to guide our expectations for the 
production of ID states in B decays, so we must rely for 
the moment on theoretical calculations [34] that suggest 
production rates roughly comparable to those for other 
charmonium states. 



TABLE III: Calculated and observed rates for radiative tran- 
sitions among charmonium levels in the potential (1). 



Transition 


k (MeV) 


Partial width (keV) 
Computed Measured" 


, Ml 

W — > '?c7 


115 


1.92 


1.13±0.41 


XcO ^ J/ip7 


303 


120 {IQ^f 


98±43 


Xci J/i/'7 


390 


242 (215)'' 


240±51 


Xc2 — y J/-07 


429 


OIK /oon\ ^ 
olo \Ia\)) 


OTn_i_ A a 


he — > Vol 


504 


482 




1 El , 
rjc — > hc'r 


126 


51 




1 El 

W — * Xc27 


128 


29 (25)' 


22ib5 


, El 

— * Xci7 


171 


41 (31) 


24zb5 


, El 

V' — >Xc07 


261 






, Ml , 

V' — >'7c7 


32 


0.04 




/ Ml 

V — »■ rtcl 


638 


n Qi 
u.y i 


n 7c:_i_n oc; 


^^(3770) ^ Xc27 


208 


3.7 




V(3770) ^ Xci7 


250 


94 




i/'(3770) ^ Xco7 


338 


287 




r,c2 ^ V (3770)7 


45 


0.34 




??c2 ^ /lc7 


278 


303 




El 

02 > Xc27 


250 


56 




El 

■02 > Xcl7 


292 


260 





"Derived from Rcf. [22]. 

''Corrected for coupling to decay channels as in Ref. [16]. 



Observing the missing narrow states. Radiative tran- 
sitions among charmonium levels are the key to dis- 
covering the remaining narrow states. Approximately 
90K B — i- Kric events are produced in the Belle ex- 
periment's data sample. Using the production rates of 
Table IV, and making the plausible assumption that 
B{B + {cc))/B{B ^ X + (cc)) is univer- 

sal, we estimate that 70K Krjc events are directly pro- 
duced, 8. IK events arise in the cascade B Kr]c2 — > 
ii:7(280 MeV)/ic ^ Js:7(280 MeV)7(500 MeV)r;c, and 
11. 7K events arise from B -» Khc 1^7(500 MeV)r]c, 
and that the sample that yielded the 39 ± 11 rj'^ discov- 
ery events was about 30K events. Likewise, the large 
radiative branching ratios of ^2 to Xci, 2 and of Xci, 2 
to J/0 provide another striking double-gamma transi- 
tion with B{B X^2) Ej=i 2 Bi^2 ^ XcJ7)S(Xc7 ^ 
J/i;-f)/B{B XJ/^P) ~ 0.12. 'The signal B Kr]c2 ^ 
i4r7(280 MeV)hc K"/{280 MeV) + hadrons may also 
provide a simultaneous observation of r]c2 and he- 

We close with a few examples of the insights to be ex- 
pected from the discovery and investigation of the miss- 
ing charmonium levels. The displacement of the l^Pi he 
from the l^Pj centroid is sensitive to Lorentz structure 
of the interquark potential. The ip'~Vc splitting is sen- 
sitive to a number of influences beyond simple potential 
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TABLE IV: Measured and estimated branching fractions for 
B decays to quarkonium levels. 



cc state r{B (cc) + X)/r{B all) (%) 

T% rfc « 0.53°' 

l^Si J/!/> 0.789 ± 0.010 ± 0.034*'= 

i^Pi he 0.132 ± o.oeo'' 

l^Po Xco 0.029 ± 0.012'* 

l^Pi Xci 0.353 ± 0.034 ± 0.024*"= 

r'P2 Xc2 0.137 ±0.058 ±0.012'' 

2^80 v'c ~ 0.18" 

2^81 V' 0.275 ± 0.020 ± 0.029*" 

l^Da Vc2 0.23-^ 

l^Di ip 0.28^ 

1^D2 %b2 0.46^ 

l"'D:l ( :i I). 65' 



"Scaled from '^Si rate. 

'Data from [32] and [33]. 

'=Known feed-down from 2S state removed. 

'^Scaled from ^Pi,2 rates using Eqn. (6). 

'Known feed-down from 2S and IP states removed. 

■'^Computed; see [34]. 

models, including the effect of virtual decay channels [37]. 
The positions of r]c2 and tp2 will further constrain ana- 
lytic calculations of spin-dcpcndcnt forces. When com- 
pared with ■!/;(3770), they will provide another test of the 
influence of decay channels in the charm threshold re- 
gion. Observation of V^(3770) in B ^ K'^*'>DD will help 
to calibrate expectations for the production of the nar- 
row states. The same final state might yield evidence for 
tp3 ^Ds. The details of V'(3770) decays are sensitive to 
S-D mixing [27]. 

Outlook: The discovery of ij'^ as a product of B decays 
realizes a long-held hope and raises new possibilities for 
filling out the charmonium spectrum. The CP violation 
experiments will enrich our knowledge of B — > (cc) + X, 
aiding our ability to estimate the production of unknown 
states. 
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